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A conventional Anger camera is not adapted to bremsstrahlung imaging and, as a result, 
even using a reduced energy acquisition window, geometric x-rays represent <15% of the 
recorded events. This increases noise, limits the contrast, and reduces the quantification 
accuracy. Monte Carlo (MC) simulations of energy spectra showed that a camera based 
on a 30-mm-thick BGO crystal and equipped with a high energy pinhole collimator is well- 
adapted to bremsstrahlung imaging. The total scatter contamination is reduced by a factor 
10 versus a conventional Nal camera equipped with a high energy parallel hole collimator 
enabling acquisition using an extended energy window ranging from 50 to 350 keV. By 
using the recorded event energy in the reconstruction method, shorter acquisition time 
and reduced orbit range will be usable allowing the design of a simplified mobile gantry. 
This is more convenient for use in a busy catheterization room. After injecting a safe activity, 
a fast single photon emission computed tomography could be performed without moving 
the catheter tip in order to assess the liver dosimetry and estimate the additional safe 
activity that could still be injected. Further long running time MC simulations of realistic 
acquisitions will allow assessing the quantification capability of such system. Simultane- 
ously, a dedicated bremsstrahlung prototype camera reusing PMT-BGO blocks coming 
from a retired PET system is currently under design for further evaluation. 
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INTRODUCTION 

Recent studies proved the tumor response and the hepatic toxicity 
correlated to their respective absorbed doses in liver 90 Y radioem- 
bolization (1-5). This strongly supports the need to individually 
determine the maximum safe activity that can be injected in order 
to optimize the patient outcome. Unfortunately, in the mean 
time, accumulating evidence showed that pre-therapy 99m Tc- 
macroagregated (MAA) distribution is not reliable in predicting 
the 90 Y microsphere distribution on a patient by patient basis 
(6, 7). Thus, a way to quantitatively assess the 90 Y microsphere 
distribution in the catheterization room is highly valuable. 

Quantitatively imaging 90 Y with an Anger camera is one of the 
most challenging topics in nuclear medicine. The bremsstrahlung 
x-rays are spread along a continuous spectrum extending up to 
the maximal beta energy emission, i.e., 2.3 MeV. The maximum 
energy usable by a camera using a mechanical collimator being 
about 0.5 MeV, acquisitions are contaminated by high energy x- 
rays scattered down into the acquisition energy window. The major 
problematic effects are: (1) the scattering inside the patient body, 
(2) the collimator scattering and penetration, (3) the lead fluo- 
rescence around 80 keV in the collimator, (4) the back- scattering 
from the light guide, the photo-multiplier tubes (PMTs) , electronic 
boards, and lead housing of the camera. As a result, contrary to 
gamma emission imaging where about 60% of the recorded events 
are geometric y-rays, geometric x-rays in bremsstrahlung imaging 
represent <15% of the recorded events (8). Numerous methods 



have been proposed to address these issues [see (9, 10) for an 
extensive bibliography] . 

State of the art bremsstrahlung single photon emission com- 
puted tomography (SPECT) implements collimator-detector 
table point spread function (PSF) pre -calculated by Monte Carlo 
(MC) in the iterative reconstruction process (11, 12). These meth- 
ods show improved quantification in phantoms (12) and in patient 
studies as well (13). However, in nuclear medicine imaging it is 
always profitable to improve the hardware performance in order 
to select the right events, rather than to correct for contaminating 
events afterward. Indeed, this last solution inevitably results in a 
higher noise level regarding the statistical nature of geometric and 
contaminating x-rays. Other choices than parallel hole collimators 
have been recently considered. 

Van Holen et al. (14) proposed the use of a rotating slat colli- 
mator that has a much higher geometric efficiency than a parallel 
hole collimator. The relative importance of septal penetration is 
reduced, resulting in a better contrast to noise ratio. Note that, 
regarding only the geometric x-rays, the high geometric efficiency 
improvement of the rotating slat collimator is counterbalanced 
by less information provided about the x-rays coming direction. 
Also, the high energy x-rays fraction crossing the crystal, and after- 
ward back scattered down in the energy acquisition window, is not 
reduced. 

Walrand et al. (15) has used a medium energy pin- 
hole (MEPH) collimator in bremsstrahlung SPECT. Collimator 
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FIGURE 1 | Gate modeling of the GE400AC equipped with the MEPH 
collimator. From right to left: blue: perpex surrounding the 90 Y point source, 
gray: 6mm-aperture tungsten insert, red: lead housing, yellow: 1/2"-thick 
Nal crystal, white: 12 mm-thick pyrex light guide, gray: 0.2 mm-thick 
hexagonal mu-metal shielding, purple: cathode-dynode cascade (PMTs 
pyrex housing was set no visible), gray: 1 cm-thick perpex modeling the 
electronics. 



penetration-scattering can only occur on the small area of the 
tungsten nose insert in contrast with parallel hole collimator 
where the scattering occurs on the inner walls of each hole. Eval- 
uation of a cold and hot spheres phantom showed that MEPH 
SPECT provided quantification accuracy similar to that of time 
of flight-positron emission tomography (TOF-PET), but with 
significantly less noise. Helical MEPH SPECTs of a liver radioem- 
bolization phantom were also acquired and showed that repro- 
ducibly accurate activity quantification can be obtained within a 
3 min acquisition time. 

In this paper, we evaluate by MC simulations the optimal design 
of an Anger camera-collimator system for bremsstrahlung imag- 
ing. We discuss the reasons why the use of the recorded event 
energy allowed by such optimal system should allow shorten- 
ing the acquisition time and reducing the camera orbit range in 
SPECT. These improvements are well- adapted to fast liver dosime- 
try assessment in the catheterization room. Last, we introduce our 
current prototype project evaluation. 

MC CODE 

Simulations were performed in the vGate 2.1 environment (16) 
running Gate 6.2 and Geant4.9.5.p01 (17). Geant4 (GEometry 
ANd Tracking) is a toolkit for the simulation of the passage of 
particles through matter initially developed in order to model 
particle and nuclear physics experiments. GATE (Geant4 Appli- 
cation for Tomographic Emission) is an advanced opensource 
software developed by the international OpenGATE collaboration 
and dedicated to numerical simulations in medical imaging and 
radiotherapy using Geant4. In order to profit from recent cor- 
rections, the source files handling the atom de- excitations and 
fluorescence emissions of Geant4.9.5.p01 were replaced by those 
of Geant4.9.6.p02. Accordingly, the call to these processes was 
updated in the LiverMore model source file. Afterward, Geant4 and 
Gate were successively rebuilt. Simulations were performed using 
the LiverMore model for the photoelectric and for the Compton 
processes. The relative yields of the Pb fluorescence peaks checked 
in the simulations were in agreement with those reported in the 
G4EMLOW (Geant4 ElectroMagnetic Low energy) data file. 

MC VALIDATION 

In order to validate our MC simulations, a 90 Y point source sur- 
rounded by 1 cm perpex was acquired using a l/2 /r -thick Nal 
(Sodium Iodine) GE400AC camera (General Electric, Milwau- 
kee, WI, USA): our only camera having a MEPH collimator. The 
GE400AC was successively equipped with a low energy high res- 
olution (LEHR), a medium energy general purpose (MEGP), a 
high energy general purpose (HEGP), and an MEPH collima- 
tor. The distance between the source and the collimators aperture 
was 10 cm. 

A retired PMT from the GE400AC was taken apart for mea- 
surement of the housing pyrex thickness, of the cathode-dynodes 
cascade dimensions and weight. The PMTs were modeled by 
hexagonal-cylindrical pyrex tubes, each one containing a cop- 
per medium of volume and of effective density derived from the 
previous measurements (Figure 1). The hexagonal bases of the 
PMTs were surrounded by 0.2 mm-thick mu-metal. A 12 mm- 
thick pyrex light guide was set between the PMTs and the crystal. 



Last, a 1 cm-thick perpex plate was set behind the PMTs in order 
to model the electronics. 

Figure 2 shows the comparison between the acquisitions and 
the MC simulations. Overall, the total counts shapes were well- 
reproduced for a wide variety of parallel hole collimators, ranging 
from a low energy to a high energy parallel hole collimator, and 
including an MEPH collimator as well. The lead fluorescence peaks 
were well-reproduced for all the parallel hole collimators and also 
its absence when using the pinhole collimator. This absence arises 
from the fact that the x-rays cannot reach the inner side of the 
MEPH housing (15). In contrary, using a parallel hole collima- 
tor, x-rays hit the inner side of all the collimator holes producing 
fluorescence photons, which are able to reach the crystal. 

Some small discrepancies were observed below the lead fluo- 
rescence peaks for the MEGP and HEGP collimators, while the 
back- scatter peak was a little bit overestimated for the LEHR col- 
limator and a little bit underestimated for the other ones. One 
possible explanation is that these collimators were produced at 
different epochs and by different companies (our GE400AC has a 
26 years long history). Manufactures often use lead alloys contain- 
ing some percentage of tin (up to 1 5% and unknown to us) in order 
to improve the casting properties. Last, for obscure reasons, the 
GE400AC energy spectra displayed structured oscillations between 
100 and 220 keV. 

OPTIMAL COLLIMATOR DESIGN: MC SIMULATIONS 

There are three adverse effects, which occur in the collimator: the 
lead K-shell fluorescence, the penetration, and the scattering (15). 

Lead fluorescence in a pinhole collimator can be reduced by 
using a bi- conical insert made in tungsten for the cone facing 
the activity, but made in tin for the cone facing the crystal in 
order to attenuate the tungsten and lead fluorescence x-rays (see 
second last section). Tin and antinomy are the highest atomic 



Frontiers in Oncology | Cancer Imaging and Diagnosis 



June 2014 | Volume 4 | Article 149 | 2 



Walrand et al. 



Camera design for bremsstrahlung imaging 




keV keV 



FIGURE 2 | Comparison between the acquisitions (black curves) and the 
MC simulations (red curves) for a 90 Y point source surrounded by 1 cm of 
perpex and set 10 cm far from a LEHR, an MEGP, a HEGP, and an MEPH 
collimator, all mounted on a 1/2" -thick Nal GE400AC camera, f is the 



septal thickness, d the hole diameter, / the hole length. For MEPH, T is the 
front housing thickness, and f the focal length. Note the absence of lead 
fluorescence peak for the MEPH collimator. The oscillations observed in the 
acquisitions from 100 to 220 keV (17 keV-period) were not explained. 



number elements (Z = 50 and 51, respectively) with high den- 
sity (~7kg/dm 3 ) having their K-shell fluorescence below 50keV 
[~30 keV (18)]. Additionally, the inner side of the lead front plate 
of the pinhole collimator can be covered by a 2mm-thick tin 
layer in order to attenuate the fluorescence produced by the small 
amount of high energy x-rays crossing the lead plate. Except for 
the GE400AC, all the pinhole collimators were modeled with these 
two features. 

Figure 3 shows the penetration and scatter component of var- 
ious collimators, i.e., the x-rays passing through the collimator 
material without any interaction or with Compton scattering. 
Most of these x-rays are above the usual energy acquisition win- 
dow. However, some of them will leave some energy within the 
acquisition window by Compton scattering in the crystal and some 
will be back scattered down into the energy acquisition window 
by the camera compartment. This dramatically hampers the image 
contrast and the quantification accuracy (9). 

The penetration and scatter component is reduced when using 
a HEGP collimator in place of the MEGP one. However for a 
similar spatial resolution, a conventional MEPH collimator (2 cm- 
housing wall thickness) already displays a lower penetration and 
scatter component. Furthermore, the septa thickness of a parallel 
hole collimator cannot indefinitely be increased without hamper- 
ing the spatial resolution and sensitivity. In contrary, the housing 
wall thickness of a pinhole collimator can be increased without 
any drawback on the collimator performances, the holding weight 
capacity of the gantry being the only limitation. A 45 cm x 25 cm 



front lead housing wall, a size well-adapted to whole-liver imag- 
ing, weights about 64 kg for a 5 cm-thickness. Further increasing 
of the thickness only resulted in small improvement (data not 
shown). 

OPTIMAL CRYSTAL CHOICE 

The main drawback of a conventional Anger camera for 
bremsstrahlung imaging is the back- scattering by the pyrex light 
guide and PMTs of the high energy x-rays passing through the 
crystal (9). 

Heard et al. (8) assessed by MC simulation the lead fluores- 
cence, the geometric x-rays, and the camera back-scatter x-rays 
originating from a 90 Y point source in water detected by a 3/8"- 
thick Nal Anger camera equipped with a medium energy parallel 
hole collimator. Below 90 keV, due to the limited crystal energy 
resolution, the geometric x-rays are significantly contaminated by 
lead fluorescence and above 90 keV the geometric x-rays repre- 
sent only a decreasing fraction of the camera back-scatter x-rays. 
As a result, the energy window optimizing the contrast is often 
constrained to (100, 150) keV (i.e., between 100 and 150 keV) (8), 
which limits the sensitivity. 

The camera back- scatter contamination can be reduced by 
using a thick crystal having a high photoelectric efficiency. This 
solution also has the additional benefit of improving the geomet- 
ric x-rays sensitivity. Gd2Si05 (GSO, gadolinium oxyorthosilicate) 
and Bi4Ge30i2 (BGO, bismuth germanate) are two good candi- 
dates. BGO has an absorption efficiency a little better than that 
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FIGURE 3 | MC simulations of the x-rays counts passing through the 
collimator material (blue curves) for a 90 Y point source surrounded by 
1 cm perpex and set 10 cm far from the collimator. The curves were 
rescaled in order to get the same geometric x-rays component (black curve) 
for all the collimators. Septa thickness is 1.4 and 3.2 mm, hole diameter is 
3.4 and 4.0 mm, hole length is 42 and 40 mm, for the MEGP and HEGP 
collimator, respectively. The length in centimeters is the thickness of the 
lead front plate of the pinhole (PH) collimator. The aperture diameter is 
8 mm and the focal length is 15 cm. 



and 320 keV using Nal and BGO, respectively. Using a 5 cm-thick 
pinhole collimator, this last intersection was shifted to 450 keV. 

POTENTIAL BENEFIT PROVIDED BY THE EXTENDED 
ACQUISITION ENERGY WINDOW 

In addition to improve the sensitivity, the extension of the acqui- 
sition energy window could provide a much more important 
benefit. 

Let's assume that we have an ideal camera-collimator system, 
i.e., no fluorescence, no scattering, no septal or housing wall pene- 
tration, a perfect energy, and a perfect spatial resolution, and that 
the patient tissue in the liver slice is similar to water (which is 
a reasonable assumption in liver imaging). Then, neglecting the 
scattered x-rays, the count rate c(x, y, E) recorded at the energy E 
in the pixel (x, y) of the camera stationary above the patient is: 

poo 

c(x,y,E) = / a(E)A(x , (l) y y , (l) y z , (l))G(x , (l) y y , (l) y z , (l)) 
Jo 

e -\LW(!-ip(*>y)) dl (l) 

where A(x f , y' , z') is the 90 Y activity, G(x\ y' z') the geomet- 
ric efficiency of the collimator, a(E) is the bremsstrahlung x-rays 
yields in water at energy E, \x(E) the x-ray attenuation coefficient 
of water at energy E. [x f (l), y f (l), z'(l)} is the parametric equation 
of the collimator line of response (Figure 6), i.e., for a parallel hole 
collimator: 



of GSO: 90 and 84% for 30 mm -thickness at 511 keV, respec- 
tively, both being much better than that of a conventional 3/8 /r - 
thick Nal crystal, i.e., ~15% (Figure 4). GSO has the benefit to 
produce a little bit more light than BGO, i.e., 12.5 and 9 pho- 
tons/keV, respectively, providing a better energy resolution. BGO 
has the drawback of a strong temperature dependence of the 
light production, i.e., — 1.2%/K. This can induce some shift of 
the energy spectrum. However, when the purpose is to acquire 
bremsstrahlung x-rays within a large energy acquisition win- 
dow, energy resolution and temperature dependence have only 
a marginal impact. 

The best crystal for a dedicated bremsstrahlung camera is 
CdW04 (cadmium tungstate) that has an absorption efficiency 
a little better than BGO (Figure 4) and a light yield slightly better 
than GSO, i.e., 13.5 photons/keV. However, the primary light decay 
of CdW04 is slow (14 |xs). As a result, the field of view should be 
split in independent crystal-PMT-application specific integrated 
circuit (ASIC) units in order to handle the count rate required in 
liver radioembolization. 

BGO CAMERA DESIGN: MC SIMULATIONS 
Figure 5 shows the total scatter component, i.e., collimator scatter- 
penetration + camera backscattering, using MEGP, HEGP, and 
high energy pinhole (HEPH) collimator mounted on a 3/8" -thick 
Nal and on a 30 mm-thick BGO camera, both with a crystal size of 
45 cm x 25 cm. Simulations were performed with a 6.6 cm-thick 
slab of 66% density Pyrex modeling, the camera back compart- 
ment (8). The 90 Y point source surrounded by 1 cm of perpex 
was set 10 cm far from the collimator. With a HEGP collimator, 
the scatter component crossed the geometric x-rays around 150 



X \l) =X .y\l )= y. Z \l) = l 

and for a pinhole collimator: 



(2) 



*'(/) 



z'{l) = F + 



y?F 2 + x 2 +y 2 " / K ' y/F 1 + x 2 + y 2 
IF 



,/F 2 + x 2 +y 2 



(3) 



/ is the distance on the line of response from the crystal for the 
parallel hole collimator and from the aperture for the pinhole col- 
limator. lp(x, y) is the / value on the intersection of the line of 
response with the patient surface facing the collimator. 
Dividing both sides of Eq. 1 by a(E)e lx(E)l P (x ' y \ we get: 

C(x,y, E) 
poo 

= / A(x'(l),y'(l),z'(l)) G(x\l),y\l),z\l))e-^ E)x dl 
Jo 

(4) 



where 



C(x,y y E) = c(x y y,E)e- [l{E)l P {x > y) /a(E) 



(5) 



|x(E) for water is a bijective function ranging from 0 (£ = oo) to 
oo (£ = 0) (18). Inverting this function gives: 



poo 

C(x,/,|x) = / A(l)e~^ l dl 
Jo 



(6) 
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FIGURE 5 | Comparison of the total scatter components (colored 
curves) with the geometric x-rays (black curves) obtained by MC 
simulations for a 90 Y point source surrounded by 1 cm of perpex and 
set 10 cm in front of the collimator. Septa thickness is 1.4 and 3.2 mm, 

hole diameter is 3.4 and 4.0 mm, hole length is 42 and 40 mm, for the 
MEGPand HEGP collimator, respectively. Front wall thickness of 5 cm, 
aperture diameter of 8 mm, focal length of 15cm for the HEPH collimator. 
45cm x 25cm crystal: 3/8"-thick Nal and 30mm-thick BGO.The purple 
curve corresponds to the BGO prototype (see second last section). The 
BGO to Nal geometric ratio is 4.8 at 400 keV, which is in line with the 4.9 
predicted by previous MC simulations (Figure 4). 



where: 

C(x y y, |x) = C(x, y, E([i)) (7) 
A(l) = A(x\l\y\l) y z r {l))G{x\l) y y\l) y z{l)) (8) 

Equation 6 is just the Laplace transform, the inversion of which 
arises in various engineering and physics problems, e.g., such as 
electronic circuits analysis or heat conduction modeling (21). The 
inversion is unique except for an additional arbitrary null-function 
N(Z), i.e., fl N(l')dV = 0VJ > 0, which is known as Lerch's the- 
orem. For the engineers and physicists, N(l) may almost be taken 
as zero (21). Thus, in theory, Eq. 6 allows us to compute A(x f , 
y\ z') via Eq. 8. Equation 6 is the equivalent of the Radon trans- 
form where the energy £, via the attenuation, plays the role of the 
rotation angle cp in conventional SPECT. 

However, the inversion of the Laplace transform appears to be 
a much more ill posed problem than conventional SPECT (22), 
especially when the Laplace transform is known, or is measured, 
only on the real positive axis (23, 24). Numerous inversion meth- 
ods and regularization schemes exist, and the choice of the optimal 
couple requires a detailed analysis of the problem to be solved, 
supported by numerous simulations. 

In real situation, similarly to conventional SPECT, the acquisi- 
tion will be affected by the limited energy and spatial resolution 
of the camera, and by the Compton scattering occurring inside 
the camera-collimator and inside the patient as well, which can 
be simulated only by MC methods. Simulation of a single acquisi- 
tion requires some months of CPU time on a 20 cores computer. 
A faster MC algorithm dedicated to bremsstrahlung is under 
development. 

Likely, regarding the noise inherent to radioactive decay imag- 
ing, Eq. 6 alone will not be sufficient to obtain a reliable activity 
reconstruction. However, similarly to the time of flight (TOF) 
introduction in PET (25, 26), the utilization of the recorded event 
energy in the reconstruction process should allow a reduction of 



the acquisition time and on the angular orbit range, resulting in a 
simplified camera mobile gantry. This is important for the use in 
a busy catheterization room. After injecting a safe activity, a fast 
SPECT, without moving the catheter tip, could be performed in 
order to assess the liver dosimetry and estimate the additional safe 
activity that could still be injected. 
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FIGURE 6 | Coordinates convention for the pinhole collimator line of 
response. 



BGO BREMSSTRAHLUNG DEDICATED CAMERA PROTOTYPE 

In order to assess the performance of a dedicated bremsstrahlung 
camera-collimator and of the benefit of the extended acquisition 
energy window, we plan to build a camera prototype by reusing 
PMT-BGO blocks of a retired Exact HR + PET (CTI, Knoxville, 
TN, USA) . The BGO block thickness is 30 mm with an overall front 
size of 35 mm x 38 mm divided into 8x8 pixels. The prototype 
will be made of 1 1 rows, each having 7 PMT-BGO blocks set along 
an arc (Figure 7) in order to avoid gaps between the front sides of 



the BGO blocks that have a trapezoidal shape. The detector FOV 
will be 25 cm x 43 cm. The detector arc shape is useful when using 
the camera with a pinhole by reducing the parallax error in the 
rows direction due to the depth of interaction. 

Similarly that for the GE400AC modeling, a BGO -block was 
disassembled and one PMT was taken apart to measure the pyrex 
and cathode-dynode cascade dimensions and weight. An accurate 
modeling of the prototype equipped with a pinhole collimator 
(5 cm- thick lead front plate, 8 mm-aperture tungsten insert) was 
performed (Figure 7). The MC simulations showed that an energy 
acquisition window extended to (50, 350) keV could be used 
(Figure 5). Access to a super computer will be required in order to 
simulate a complete acquisition. 

In order to account for the different interactions of the x-rays 
in the collimator and camera, the PSF (x, y, z, E) of a 90 Y point 
source will be measured for different water depths z, and used 
in the reconstruction process. Since the BGO light yield is tem- 
perature dependent (— 1.2%/K), use of the recorded event energy 
with BGO requires measuring these PSF with the camera being 
in thermal equilibrium. The energy spectrum of each block for a 
99m Tc source will also be recorded in this thermal equilibrium. A 
thin temperature probe will be set in contact with the two central 
BGO blocks to assess the thermal equilibrium. Before each patient 
acquisition, the energy spectrum of the blocks using a 99m Tc source 
will again be recorded and used to rescale the energy measurement 
of the bremsstrahlung x-rays during the liver acquisition in order 
to account for difference in temperatures from 1 day to another 
one. 

CONCLUSION 

Monte Carlo simulations prove a significant reduction of the 
total scatter component when using a 30mm-thick BGO cam- 
era equipped with a HEPH collimator (i.e., 5 cm-thick). With this 
camera, the total scatter component cross the geometric x-rays 
counts around 450 keV rather than below 150 keV when using 
a conventional 3/8 /r -NaI camera equipped with a parallel hole 
collimator. 

In addition to improve conventional SPECT, mathematics 
shows that such novel performances should allow the develop- 
ment of a simplified mobile gantry, which is more convenient for 
use in a busy catheterization room during the liver radioemboliza- 
tion. Liver dosimetry assessment will be possible without moving 
the catheter tip in order to estimate the additional safe activity that 
could still be injected. The development of a faster MC code dedi- 
cated to bremsstrahlung imaging is under way in order to predict 
the feasible camera orbit range reduction. 
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FIGURE 7 | Gate modeling of the prototype made of 7 x 11 BGO-blocks 
coming from a retired Exact HR + PET and equipped with a pinhole 
collimator (5 cm-thick lead front plate, 8 mm-aperture diameter). From 
right to left: blue: 1 cm-thick perpex surrounding the 90 Y point source, red: 
lead housing, white: tungsten outer conical insert, blue: tin inner conical 
insert and 2mm-thick layer, yellow: 30mm-thick BG, gray: mu-metal, 
purple: cathode-dynodes cascade (PMT pyrex housing walls were set no 
visible), black: PMT plastic sockets, white: perpex modeling the electronics. 
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